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SECTION  1 


INTRODUCTION 


During  recent  years,  most  of  the  studies  of  high  altitude  nuclear 
bursts  have  emphasized  the  effect  of  the  nuclear  environment  on  radar  and 
optical  components  of  weapon  systems.  Chemistry  calculations  done  as  a 
part  of  these  studies  generally  assumed  a substantial  initial  disturbance 
of  ambient  conditions.  It  was  also  generally  assumed  that  the  nuclear 
effects  were  "over"  when  the  electron  density  had  decayed  to  approximately 

106  cm"3.  Natural  ionization  sources  which  are  negligible  in  comparison 

« 

with  the  nuclear  sources  were  riot  considered. 

There  is  now  interest  in  extending  our  current  capability  to 
compute  high  altitude  nuclear  effects  to  include  results  that  can  be  used 
in  communication  system  studies.  One  of  the  many  questions  that  needs  to 
be  addressed  concerns  the  adequacy  of  the  deionization  chemistry  calculations. 
Whereas  it  is  usually  assumed  that  radars  would  operate  satisfactorily  in 
a non-nuclear  environment,  it  is  well-known  that  communication  links  over 
a large  part  of  the  frequency  spectrum  can  be  substantially  enhanced  or 
degraded  by  naturally  occurring  variations  in  ionospheric  conditions.  We 
therefore  need  to  verify  that  our  calculations  of  the  ionization  and  de- 
ionization that  occur  following  a nuclear  burst  limit  correctly  to  the 
ambient  electron  and  ion  density  in  the  absence  of  a nuclear  source.  This 
in  turn  implies  that  natural  ionization  sources  can  no  longer  be  ignored. 

This  report  describes  our  efforts  to  incorporate  natural  ionization 
sources  into  the  chemistry  codes  that  are  used  to  compute  the  deionization 
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following  a nuclear  burst.  We  then  compute  the  density  of  electrons  and 
other  species  of  interest  in  an  ambient  atmosphere  and  compare  the  results 
with  observations. 
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SECTION  2 


NATURAL  IONIZATION  SOURCES 


DAYTIME  SOURCES 

The  primary  ion  source  in  the  daytime  ambient  ionosphere  is  solar 
radiation  in  the  ultraviolet,  extreme  ultraviolet  and  x-ray  regions  of  the 
spectrum.  Molecular  oxygen  can  be  ionized  by  radiation  with  wavelengths 
less  than  1028  k,  while  the  molecular  nitrogen  ionization  threshold  occurs 
at  796  A.  The  atomic  oxygen  and  nitrogen  thresholds  are  at  910  A and  853  A 
respectively.  Since  atomic  nitrogen  is  a very  minor  species  at  all  altitudes, 
its  photoionization  can  be  safely  neglected  under  ambient  conditions  and 
none  of  the  calculations  reported  here  include  nitrogen  atom  photoionization 

Another  minor  species,  NO,  cannot  be  neglected  as  its  ionization 
threshold  occurs  at  1338  A and  a strong  peak  in  the  solar  flux  associated 
with  hydrogen  Lyman-a  (H  Ly-a)  emission  occurs  at  1215  A.  This  H Ly-a 
radiation  is  not  strongly  absorbed  by  the  major  atmospheric  species  and 
penetrates  to  lower  altitudes  than  the  radiation  that  ionizes  N2,  02,  and 
0.  In  the  D-region  N0+  is  the  primary  ion  produced  by  absorption  of  solar 
radiation.  A discussion  of  ionization  produced  at  lower  altitudes  by  solar 
radiation  has  been  given  by  Bourdreau,  et  al1. 

The  calculation  of  a solar  photoionization  rate  constant  at  an 
altitude,  z requires  a knowledge  of 

1.  The  solar  flux  at  the  top  of  the  atmosphere, 
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2.  absorption  and  ionization  cross-sections  for  the  important 
atmospheric  species, 

3.  densities  of  the  atmospheric  species  along  the  absorption 
path  of  the  solar  ray  from  the  topside  down  to  the  altitude  z. 

Such  a calculated  ionization  rate  constant  will  thus  reflect,  in  varying 
degrees,  the  level  of  solar  activity  which  in  turn  is  linked  to  atmospheric 
variations  in  temperature,  and  density  associated  with  latitude,  season,  etc. 

Detailed  topside  measurements2  of  the  solar  spectrum  are  now 
available  from  satellite  data,  and  tabulations  of  cross-sections  and  standard 
atmospheric  density  profiles  can  be  found  in  a variety  of  sources3*5  . 

Several  computations  of  solar  photoionization  rate  constants  have 
also  been  published1*-6  . For  the  work  reported  here,  we  have  adopted  the  set 
of  rate  constants  given  by  Keneshea,  et  al.,  in  Chapter  13  of  Reference  5 
for  N2,  02,  0 and  NO  as  function  of  altitude  and  solar  zenith  angle. 

These  rate  constants  represent  conditions  of  low  solar  activity  and  correspond 
to  a latitude  of  30°.  Other  sets  of  rate  constants  appropriate  to  different 
conditions  could  easily  be  added  by  expanding  the  tabular  input  to  the 
chemistry  codes. 

Cosmic  rays  penetrate  deeply  into  the  atmosphere  causing  ionization 
proportional  to  the  atmospheric  density.  Thus,  above  60-70  km,  cosmic  rays 
as  a daytime  ion  source  can  generally  be  ignored.  Furthermore,  the  intensity 
of  the  cosmic  rays  decreases  with  increasing  solar  activity  and  decreasing 
latitude.  The  existing  D-region  chemistry  code  contains  a background  cosmic 
ray  source7  and  this  has  been  retained.  The  E and  F-region  codes  do  not 
contain  such  a source  and  none  was  added. 
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NIGHTTIME  SOURCES 


In  the  E and  Fi  regions  nighttime  ionization  levels  can  be  several 
orders  of  magnitude  below  daytime  levels.  Such  levels  cannot  be  explained 
by  simply  removing  the  daytime  source.  Swider's8study  of  nighttime  ionospheric 
processes  indicates  that  some  source  is  required  to  maintain  nighttime  E 
region  levels  and  he  suggested  that  possible  sources  might  be  scattered 
hydrogen  Ly-a  and  Ly-8  radiation.  Recent  rocket  probes  into  the  upper 
atmosphere  have  confirmed  the  existence  of  H Ly-a  and  H Ly-8  nightglow 
as  well  as  nightglow  due  to  He (504  A)  and  HE+(389  A)  lines9.  The  resulting 
data  indicate  that  there  is  sufficient  flux  of  scattered  radiation  to  account 
for  the  observed  maintenance  of  the  nighttime  ionosphere10. 

Following  Keneshea,  et  al.11,  we  have  attempted  to  simulate  a 
nighttime  source  in  the  E and  F regions  by  reducing  daytime  photoionization 
rate  constants  in  such  a way  as  to  roughly  reflect  observed  nighttime  fluxes. 

D-region  ionizing  sources  at  night  are  primarily  cosmic  rays. 

Other  possible  sources  at  night  are  the  diffuse  x-ray  background,  galactic 
sources  of  x-rays12  and  precipitating  energetic  electrons13.  Scattered 
H Ly-a  may  also  play  a role12.  Other  than  cosmic  rays,  none  of  these 
sources  has  been  included  in  the  D-region  code  at  present,  and  no  nighttime 
D-region  calculations  were  done. 

OTHER  SOURCES 

Other  possible  natural  ionization  sources  include  energetic 
electron  showers  and  proton  and  a-particle  bombardments.  These  sources  are 
generally  more  important  at  high  latitudes  and  are  often  associated  with 
disturbed  ion  levels  such  as  are  observed  in  auroral  and  other  related 
phenomena.  At  present  no  provision  for  these  types  of  sources  has  been 
made  in  our  nuclear  effects  chemistry  codes.  Chemistry  codes  oriented 
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towards  the  study  of  auroral  data  have  been  developed  at  MRC  in  support  of 
the  ICECAP  program14.  Presumably  some  of  the  techniques  and  methods  used 
in  these  studies  can  be  applied  to  our  existing  codes  should  the  inclusion 
of  natural  particle  bombardment  sources  be  required  in  future  work. 
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SECTION  3 

CHEMISTRY  CODES 


THE  D- REG I ON  (60-100  km) 

A comprehensive  code,  DCHEM,  containing  more  than  40  species 
and  several  hundred  reactions  has  been  developed  by  Scheibe15.  Included 
in  these  reactions  are  photodissociation  and  detachment  processes.  Addition- 
ally the  code  normally  expects  as  input  some  initial  high  level  ionization 
source.  For  the  calculations  reported  here,  this  source  has  been  set  to 
zero.  The  background  cosmic  ray  source  is  retained  in  its  present  form 
of  p = p*10~ 1 7 ion  pairs/cm3sec,  where  p is  the  production  rate  and 
p is  the  density.  The  solar  photoionization  rate  constants  of  Reference  5 
were  added  to  the  program. 

THE  E AND  F REGIONS  (100-300  km) 

Two  chemistry  codes  are  available  for  these  altitudes  and  they 
reflect  different  levels  of  sophistication  and  speed.  SDCCHEM  is  a relatively 
simple  code  containing  12  reactions  and  8 species,  N2,  02,  NO,  N,  0,  N+,  0+ 
and  X*  , which  represents  the  sum  of  the  molecular  ions,  NO+,  0*  and  N*. 

A full  description  of  this  code  is  contained  in  Reference  16.  The  second 
code,  RK3,  contains  15  species  including  excited  N(2D)  and  0+(2D)  and  it 
has  about  50  reactions.  Further  details  on  this  code  are  given  in  Appendix  A. 

Each  of  these  codes  contains  provisions  for  instantaneous  high 
level  ionizing  sources.  These  were  set  to  zero  for  all  the  ambient  calculations 
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reported  here  and  the  photoionization  rate  constants  of  Reference  5 were 
added. 

INITIAL  AND  FINAL  DENSITIES 

The  DCHEM  code  was  run  with  ambient  temperature  and  density  profiles, 
appropriate  to  the  equilibrium  daytime  D-region  composition  and  to  the  various 
photo-chemical  rate  constants  contained  in  the  program.  Thus,  initial 
profiles  for  C02>  H20,  0,  0^,  N02>  NO,  and  O^A)  were  included  in 
addition  to  the  major  species  densities. 

Both  E and  F region  codes  were  run  with  the  same  temperature  and 
density  profiles  used  to  compute  the  solar  photoionization  rate  constants5. 

The  inclusion  of  an  initial  NO  profile  was  found  to  be  important  in 
attaining  equilibrium.  This  profile  is  still  subject  to  a great  deal  of 
uncertainty17. 

Once  the  initial  densities  and  temperatures  were  specified,  each 
code  was  allowed  to  run  out  to  about  104  seconds.  Equilibrium  with  respect 
to  ionization  levels  was  taken  to  be  reached  whenever  the  electron  density 
remained  at  a steady  value.  At  most  altitudes  this  was  reached  between 
5*102  and  5*103  seconds.  In  some  cases  no  steady  value  was  attained.  This 
was  true  for  the  D-region  results  at  60  and  65  km  and  also  for  some  of  the 
higher  altitude  F-region  calculations. 

The  electron  density  profiles  reported  here  are  the  equilibrium 
values  as  defined  above.  If  no  equilibrium  was  reached,  the  maximum  values 
are  given  except  for  the  60  and  65  km  results.  For  these  two  altitudes 
the  computed  electron  density  begins  to  oscillate  at  late  times.  Therefore, 
the  values  reported  for  these  two  altitudes  are  the  electron  densities 
computed  at  103  seconds,  the  time  at  which  the  D-region  code  indicates  that 
equilibrium  is  first  reached  with  respect  to  the  electron  density. 
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None  of  the  electron  densities  reported  here  is  expected  to 
differ  by  more  than  a factor  of  two  or  three  from  the  true  equilibrium 
value  the  codes  would  predict  if  running  time  had  been  extended.  Experimental 
measurements  can  easily  differ  by  more  than  this.  The  experimental  profiles 
chosen  for  comparison  were  selected  to  match  the  theoretical  conditions  as 
closely  as  possible  with  respect  to  solar  activity  (low),  solar  zenith  angle 
(X*30,  60  and  75  degrees)  and  latitude  (30  degrees). 


SECTION  4 


RESULTS  AND  DISCUSSION 


THE  DAYTIME  D- REGION 

Reported  measurements  of  D-region  electron  density  profiles  show 
variations  of  as  much  as  an  order  of  magnitude  at  a given  altitude18.  These 
variations  reflect  different  levels  of  solar  activity,  season,  latitude, 
etc.,  as  well  as  different  experimental  techniques.  As  the  altitude  decreases, 
experimental  difficulties  in  detecting  the  lower  electron  densities  tend 
to  increase. 

Two  daytime  D-region  calculations  were  done  using  the  DCHEM  code. 

One  is  for  a solar  zenith  angle  of  60°  and  the  other  is  for  30°.  The 
calculated  electron  densities  are  compared  with  experimental  profiles  in 
Figures  1 and  2.  As  expected  the  computed  electron  densities  are  somewhat 
larger  when  X“30°  than  when  x=60° . 

Figure  1 shows  the  60°  results.  Two  experimental  profiles  are 
given  for  comparison.  The  Mechtly  and  Bowhill19  data  are  based  on  five 
rocket  probes  launched  from  Wallops  Island,  Virginia  (latitude  38°N)  during 
quiet  sun  years  at  X=60°.  Mitra  and  Chakrabarty's1 8 profile  is  a composite 
of  a large  number  of  mid-latitude,  low  solar-activity,  xs60°  experiments. 

Thus,  the  experimental  conditions  match  the  assumed  theoretical  conditions 
fairly  closely.  The  x=30°  profiles  in  Figure  2 are  less  well  matched,  since 
this  zenith  angle  will  be  attained  only  for  a limited  time  during  the  year 
at  mid- latitudes.  The  experimental  profile18  in  Figure  2 represents  low 
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latitude  and  medium  solar  activity  for  x - 30° . The  comparable  low-solar- 
activity,  mid- latitude,  30°  profile  would,  perhaps,  exhibit  somewhat  reduced 
electron  densities  above  65  km  and  enhanced  densities  below  this  altitude19. 


The  overall  agreement  between  both  computed  and  measured  profiles 
is  satisfactory,  considering  the  variety  of  factors  which  can  affect  measured 
electron  densities.  Mechtly  and  Bowhill19  indicate  that  their  X=60°  data 
is  accurate  to  about  10%  at  90-95  km,  but  below  these  altitudes,  the  measure- 
ments are  subject  to  increasing  uncertainties  and  error  factors  as  large 
as  two  or  three  may  occur  around  50-60  km. 

The  largest  discrepancies  between  theory  and  experiment  are  in  the 
75-85  km  region.  In  this  region  detachment  of  0~  and  0^  occurs  either 
by  collisional  association  with  0 or  by  photodetachment.  Banks4  has 
suggested  that  the  chemical  process  involving  0"  and  0~  are  more 
efficient  than  photodetachment  in  the  D-region  and  the  calculated  DCHEM 
destruction  rates  do,  in  fact,  exhibit  this  behavior.  An  excess  of  atomic 
oxygen  can  thus  increase  the  computed  electron  density.  DCHEM  consistently 
computed  late-time  atomic  oxygen  profiles  that  are  higher  than  the  initial 
concentrations  and  it  has  done  so  in  the  past  starting  from  more  disturbed 
conditions20 . 


The  most  important  source  of  0 and  0^  is  the  photodissociation 
of  0^  (the  resulting  0 atoms  react  with  0^  to  produce  the  0^) . The  rate 
constants  of  most  photodetachment  processes  occurring  in  the  D-region  are 
not  strong  functions  of  the  solar  zenith  angle,  since  significant  absorption 
of  dissociating  radiation  does  not  occur  in  or  above  the  D-region.  The 
statement  does  not  apply,  however,  to  the  0^  dissociation.  Owing  to  the 
large  density  of  02,  substantial  absorption  of  the  dissociating  radiation 
does  occur  both  in  and  above  the  D-region.  For  example  at  80  km,  the 
photodissociation  rate  for  X=60°  is  computed  to  be  about  half  of  x=0°  (overhead 
sun)  value21.  In  DCHEM,  all  photodissociation  rates  are  for  x=0°.  A 
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decrease  in  the  02  photodissociation  rate  would  decrease  atomic  oxygen 
desnity  and  the  associative  detachment  rates  for  0~  and  0"  . However, 
our  xs30°  calculation  also  tends  to  overestimate  the  late-time  atomic 
oxygen  profile,  so  that  the  inconsistencies  cannot  be  resolved  solely  on 
this  basis,  since  the  X*30°  photodissociation  rate  constant  is  essentially 
the  same  as  the  overhead  sun  value. 

The  overestimation  of  the  atomic  oxygen  may  also  be  due  to  the 
lack  of  transport  effects  in  DCHEM.  Theoretical  studies  by  Colgrave,  et  al.22 
have  indicated  that  observed  oxygen  profiles  in  the  ambient  D-region  can  be 
explained  only  by  inclusion  of  transport  effects.  Similar  conclusions  were 
reached  by  Strobel,  et  al.23,  for  NO.  A general  discussion  of  how  dynamic 
processes  can  affect  the  ionization  balance  in  the  D-region  can  be  found 
in  Reference  24. 

Most  likely  no  one  factor  can  be  isolated  to  account  for  the 
apparent  overestimation  of  the  electron  density  by  DCHEM  in  the  middle 
D-region.  We  have,  so  far,  focused  on  one  possible  source  of  "extra" 
electrons  - the  atomic  oxygen  profile.  Another  factor  is  the  uncertainty 
in  the  nitric  oxide  profile.  A lowering  of  the  NO  profile  and  hence  a 
lowering  of  the  source  of  the  primary  N0+  ions  and  electrons  by  a factor 
of  two  can  easily  be  justified25. 

Still  another  factor  is  the  question  of  ion  clustering.  Measure- 
ments made  in  the  quiescent  upper  D-region  indicate  that  both  the  positive 
and  negative  ions  are,  to  a large  extent,  hydrated;  i.e.,  they  have  formed 
clusters  with  water  molecules25"27.  The  DCHEM  code  does  not  have  negative 
ion  hydrates  and  since  these  are  probably  less  easily  detached  by  0 and 
02(3A)  the  DCHEM  code  will  tend  to  overestimate  electron  detachment.  DCHEM 
does  contain  positive  ion  clustering  but  the  ions  it  predicts  as  dominant 
are  not  those  which  are  observed.  The  explanation  for  this  is  not  yet  known. 
Those  that  are  observed  are  more  heavily  hydrated  and  have  larger  coefficients 


16 


rnmmmmmmmmmm 


•mm.  *>  ‘mrn&nmmvw  m 


of  dissociative  recombination  with  electrons.  Thus  the  DCHEM  code  under 
estimates  electron  recombination. 

All  of  these  factors  probably  combine  to  yield  the  higher  than 
observed  electron  densities  reported  here. 


THE  DAYTIME  E AND  Fx  REGIONS 

{ 

} 

Daytime  electron  density  profiles  computed  with  solar  photoionization 
rate  constants  corresponding  to  zenith  angles  of  30,  60  and  75°  are  shown 
in  Figures  3,  4,  and  5.  The  experimental  densities  at  X=406  and  60°  are 

< 

composite  mid-latitude,  low  activity  profiles  obtained  from  a large  number 
of  observations28.  The  X=75°  results  are  compared  with  rocket  data  measured 
one  hour  after  sunrise  at  Wallops  Island  in  September  1965,  a year  of  low 
solar  activity29.  All  six  calculations  reproduce  the  experimental  profiles 
quite  well  down  to  90  km. 

I 

l 

I 

As  the  figures  jhow,  the  computed  SDCCHEM  electron  densities  are 
always  lower  than  those  obtained  from  RK3. . In  the  E and  F,  regions,  electrons 
are  lost  mainly  through  dissociative  recombination  of  N0+  and  0*  . SDCCHEM 
lumps  these  two  processes  together  with  recombination  and  the  rate 

constant  for  all  three  recombination  is  controlled  by  the  fastest  processes, 
which  is  the  NO+  dissociative  recombination.  Hence  SDCCHEM  tends  to 
recombine  electrons  more  quickly  than  RK3,  which  considers  each  process 
separately. 


Figure  6 compares  the  calculated  0 and  molecular  ion  profiles 


(X2=N0  +02+N2)  at  X=60°  with  Johnson's  composite  daytime  profile  for 


low  solar  activity30.  The  computed  profiles  for  0 and  X-  reproduce 


the  relative  behavior  of  the  corresponding  experimental  densities  rather 
nicely.  Johnson's  profiles  also  include  N~  which  is  not  shown  here.  This 


ion  tends  to  remain  about  an  order  of  magnitude  less  than  the  major  X2  ions. 
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Solid  Line  - Observed,  Re 
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Electron  Density  (cm-3) 

Figure  3.  Comparison  of  Computed  and  Observed  E-  and  Fi-Region  Electron  Density 
Profiles  for  Small  Solar  Zenith  Angle. 


Solid  Line  - Observed,  Ref  28 
Dashed  Lines  - Calculated 


Figure  4.  Comparison  of  Computed  and  Observed  E-  and  F, -Region  Electron  Density 
Profiles  for  x = 608. 


Solid  Line  - Observed,  Ref  29 
Dashed  Lines  - Calculated,  x 3 


Figure  5.  Comparison  of  Computed  and  Observed  E-  and  Fi-Region  Electron  Density 
Profiles  for  Large  Solar  Zenith  Angle. 


RK3  predicted  N^  levels  are  roughly  consistent  with  the  observed  N* 
profile.  The  ion  data  in  Johnson's  composite  profile  is  normalized  to 
different  electron  density  profiles  below  and  above  130  km,  and  the  maximum 
deviation  introduced  by  this  use  of  different  profiles  is  estimated  by 
Johnson  to  be  no  more  than  25%.  The  RK3  results  are  within  25%  of  the 
experimental  densities  at  most  altitudes,  although  they  are  consistently 
lower  throughout  the  E and  F regions.  It  is  not  easy  to  determine  whether 
this  reflects  some  inconsistency  in  the  chemistry  codes  or  simply  a "normal" 
variation  of  the  experimental  profiles.  RK3  late-time  NO  densities  are 
below  the  initial  NO  profile.  This  implies  that  either  the  RK3  loss  rate 
for  NO  is  too  large  or  that  the  NO  profile  of  Reference  6 is  too  high 
for  the  conditions  under  study.  Since  the  NO  profile  is  still  subject  to 
considerable  controversy17,  it  seems  impossible  at  this  time  to  attribute 
the  depressed  late-time  NO  levels  in  RK3  solely  to  shortcomings  in  the 
theoretical  model. 

THE  NIGHTTIME  E AND  Fi  REGIGNS 

Figure  7 shows  a nighttime  electron  density  profile  computed  by 
RK3  using  a source  of  hydrogen  Ly-a  set  to  2%  of  its  daytime  X=30°  values. 
This  corresponds  to  a flux  of  H Ly-a  radiation  of  about  6><109  photons/cm2sec 
above  100  km  as  compared  to  reported  values  of  4-6><109  photons/cmzsec. 4 
The  experimental  profiles31  are  smoothed  curves  obtained  from  ionograms  at 
Boulder,  Colorado  in  1958,  a year  of  intermediate  solar  activity. 

The  agreement  between  theory  and  experiment  is  not  as  satisfactory 
as  the  daytime  results.  The  H Ly-a  produces  N0+,  but  no  other  positive 
ion,  which  is  consistent  with  observations32  that  N0+  is  the  dominant  ion. 
However,  nighttime  NoVo*  ratios  tend  to  remain  constant  at  around  0.1 
which  implies  that  some  source  is  producing  0^  10 . Scattered  H Ly-B  and 
HE  II  (304  A)  and  He  I (584  A)  fluxes  of  between  1 and  10*10®  photons/cm2sec 
have  been  reported3.  In  order  to  reduce  the  daytime  photoionization  rate 
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constant,  an  estimate  is  needed  of  the  percentage  these  particular  lines 
contribute  to  the  total  solar  spectrum  flux.  All  of  these  lines  are  fairly 
intense  compared  to  the  solar  continuum  and  a very  rough  guess  based  on  the 
flux  tabulations  of  Banks'1  would  bo  that  they  represent  around  1/10  of  the 
total  number  of  photons  below  the  0^  ionization  threshold.  The  nighttime 
fluxes  of  these  lines  are  at  least  5xl03  less  than  daytime  levels,  so  that 
a reduction  factor  of  about  5*104  or  0.05%  might  be  appropriate  and  an 
order  of  magnitude  variation  in  this  factor  could  easily  be  possible.  A 
calculation  using  0.05%  of  the  daytime  0 photoionization  rate  constant 
and  1%  of  the  NO  rate  constant  tended  to  overestimate  the  nighttime  electron 
density  as  well  as  the  NoVo*  ratio.  Since  both  experimental 3 3 and 
theoretical34  evidence  is  now  beginning  to  indicate  that  nighttime  ion 
transport  effects  in  both  the  E and  F regions  may  be  comparable  to  chemical 
production  and  loss  rates,  we  have  not  attempted  at  this  time  to  model  the 
nighttime  ionosphere  in  a more  detailed  fashion. 

Finally,  it  should  be  noted  that  the  nighttime  electron  densities 
computed  by  SDCCHEM  differed  from  the  comparable  RK3  values  in  the  same  way 
that  the  daytime  levels  did;  ..e.,  they  were  slightly  less  than  the  RK3 
values  at  all  altitudes  and  presumably  for  the  same  reasons. 

THE  DAYTIME  F2  REGION 

Figure  8 compares  the  electron  density  computed  by  RK3  with  X=30° 
to  Johnson's  composite  profile  for  low  solar  activity30.  The  computed  values 
never  reached  equilibrium  above  220  km,  an  indication,  perhaps,  that  chemical 
time  constants  are  becoming  comparable  to  ion  diffusion  times  above  this 
altitude. 


Diffusion  in  the  upper  atmosphere  may  occur  in  response  to  a 
variety  of  external  forces  such  as  gravity,  pressure  gradients,  neutral 
atmospheric  winds  and  electric  and  magnetic  fields.  The  interaction  of  the 
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Figure  8.  Comparison 


magnetosphere  and  ionosphere  is  also  important4.  Walker  and  McElroy34 
estimate  that  transport  terms  can  become  comparable  to  chemical  production 
and  loss  terms  at  220  km  and  above  during  the  day. 

Despite  the  neglect  of  diffusion  in  our  codes,  the  agreement 
between  theory  and  experiment  is  satisfactory  up  to  the  F2  peak. 
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SECTION  5 


CONCLUSIONS 


The  addition  of  solar  photoionization  processes  to  existing 
chemistry  codes  developed  for  severely  disturbed  atmospheres  associated 
with  nuclear  bursts  can  provide  reasonably  realistic  daytime  electron  density 
profiles  for  quiet  ambient  conditions  over  an  altitude  range  of  60-250  km 
(up  to  the  F2  peak) . Some  provisions  for  transport  effects  above  200  km 
is  probably  necessary  in  order  to  correctly  compute  electron  density  profiles 
above  250  km. 

Nighttime  ambient  electron  density  profiles  for  the  E and  regions 
are  less  accurately  predicted  than  daytime  profiles  and  this  may  reflect 
the  neglect  of  ion  transport  in  the  existing  codes.  Both  the  origin  and 
intensity  of  natural  nighttime  ionization  sources  are  uncertain  and  this 
must  also  be  considered  as  a contributing  factor. 

The  simplified  but  much  faster  SDCCHEM  code  can  compute  satisfactory 
E and  F region  ambient  electron  densities.  This  is  encouraging  in  terms  of 
using  the  program  for  late-time  and  large  distance  decoupled  chemistry 
calculations  in  conjunction  with  phenomenology  codes  like  MICE  where  computer 
running  time  considerations  are  important. 

The  computation  of  electron  and  ion  density  in  ambient  and  very- 
slightly-perturbed  atmospheres  is  but  one  of  many  aspects  of  the  overall 
problem  of  nuclear  effects  on  communication  systems.  Other  aspects  of  the 
problem,  such  as  calculation  of  atmospheric  motion  and  fission  debris  position 
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over  very  large  regions  and  for  very  long  times,  and  the  effect  of  the 
perturbed  environment  on  transmission  of  information,  are  being  studied 
under  other  contracts.  Our  goal  is  to  be  able  to  make  reliable  predictions 
of  communication  system  performance  in  a nuclear  environment.  We  are  unable 
to  say  at  this  time  what  accuracy  in  the  chemistry  calculations  wi’l  ultimately 
be  required  to  meet  this  goal.  However,  our  impression  from  work  reported 
here  is  that  the  chemical  aspects  of  the  problem  are  for  the  moment  at 
least,  relatively  we  11 -understood.  Less  well-in-hand  are  the  fluid  motion 
(transport)  aspects  of  the  problem,  and  the  techniques  by  which  the  fluid 
motion  and  chemistry  are  tied  together. 
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APPENDIX  A 


E AND  F REGION  RK3  CODE 

RK3  contains  nearly  50  reactions  involving  15  species.  Table  A-l 
lists  the  reactions  and  rate  constants  included  in  RK3.  The  species  followed 
are:  N2>  02>  NO,  N02,  03>  0,  N(4S),  N(2D),  N+,  0+(4S),  0+(2D) , NO+,  N*,  0* 

and  excited  N2(l).  At  present  the  N2  vibrational  temperature  is  taken  to 
be  the  same  as  the  gas  kinetic  temperature. 

Any  or  all  of  the  initial  species  densities  may  be  specified  at 
the  start  of  a run.  The  total  mass  density  and  energy  are  also  expected  as 
input.  Provision  for  conservation  of  species  is  in  the  code  and  was  used 
in  all  the  calculations  done  here. 

The  numerical  integration  of  the  15  coupled  rate  equations  is 
carried  out  by  either  one  or  both  of  two  different  methods,  depending  on 
the  stiffness  of  the  equations.  The  first  method  is  relatively  fast  and 
works  well  for  non-stiff  to  moderately  stiff  equations.  It  is  based  on  the 
procedures  developed  by  Young  and  Boris35.  The  second  and  often  slower 
scheme  carries  out  the  integration  by  the  Gear  method36.  It  is  used  whenever 
the  Young-Boris  procedure  time  step  fell  below  some  critical  value  due  to 
severe  stiffness.  For  ambient  conditions  the  Gear  method  is  almost  never 
necessary.  This  is,  however,  not  true  for  severely  disturbed  atmospheres. 

Preceding  page  blank 
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Table  A-l.  Reactions  and  Rate  Constants  In  RIG* 


REACTION 


The  rate  constants  in  this  table  are  given  in  the  form  k=aTu  e_t/l. 

In  the  a column  read  I.O(-IO)  as  l.OxlO"10. 

This  rate  constant  is  limited  to  values  of  k for  which  1.3(-12)  <k<  I.O(-IO) 
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Table  A-l  (Continued).  Reactions  and  Rate  Constants  In  RK3. 


REACTION 

a 

C 

25) 

0 + 03  - 202 

1 * 4 (-1 1 ) 

0 

2200 

26) 

O + NO  + M - N02+M 

2.9(-33) 

0 

940 

27) 

0 ♦ NO  ■*  N02  + hv 

4.11-12) 

-1.94 

0 

28) 

o + no2  - no  + o2 

1 • 6(-l 1 ) 

0 

300 

29) 

0 + N0  -"Og  + NCS) 

5.3(-12) 

0 

20200 

30) 

N+  + NO  - NO+  + Nl**S) 

4.0(-10) 

0 

0 

31) 

N+  + NO  - N0+  + N(2D) 

4. 0(-10) 

0 

0 

32) 

N^+e  - N("S)  +hv 

(As  In 

Reference  37) 

33) 

oVs)+e  -0+(2D)+e 

2.5(-8) 

0 

38504 

34) 

0+('*S)  + M -OVD)  +M 

1.0(-12) 

0 

28504 

35) 

O^-Sj+e  - 0 + hv 

(As  In 

Reference  37) 

36) 

0+(*0)  +e  - 0+(‘S)  +e 

1.0(-8) 

0 

0 

37) 

0+(2D)  +M  ■»  0+(%S)  +M 

1.0(-12) 

0 

0 

38) 

0+(2D)  + e •»  O + hv 

(As  In 

Reference  37) 

39) 

N0+  + e * { N(%S)  + | N(2D)  +0 

6.9(-6) 

-0.5 

0 

40) 

no  + Oj  - no2  + o2 

9. 5(-l 3) 

0 

1300 

41) 

O3  + hv  •*  02  + 0 

1.0(-2) 

0 

0 i 

42) 

N2  + e - NCS)  +N(2D) 

8.3(-7) 

-0.2 

0 

43) 

N2  + °2  - 02  + N2 

8.5(-10) 

-0.5 

0 

44) 

n2  + no  * no+  + n2 

3.3(-10) 

0 

0 

45) 

n2  + 0 - 1n(-s)+1n(2d)  + no 

6.0(-10) 

0 

0 

46) 

02  + e - 2 0 

l.K-5) 

-0.7 

0 

47) 

o2  + no  - o2  + no+ 

6.3(-10) 

0 

0 

48) 

n2(1)  +m  - n2  + m 

9.97(-31 ) 

4.748 

3876 

49) 

N2(l)+e  - N2  + e 

3.20(-16) 

1.5 

0 
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